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(C)  ^SITUATION,  MISSION,  AND  REPORTS  OF  RESULTS  (U) 

The  CHURCH  OPAL  Exercise  is  one  of  a series  of  LRAPF 


0. 

(C) 

Exercises  designed  specifically  to  acquire  environmental  acoustic 
data  required  for  ASW  program  decisions.  At-*sea  operations  vere 
conducted  during  September  1975  in  a region  of  the  Northeast 
Pacific  Ocean^ shown  in  Figure  0-1. 


(C)  This  exercise  includes  the  following  ocean  acoustic  ex- 

periments related  to  the  Moored  Surveillance  System  (MSS) , the 
Sound  Surveillance  System  (SOSUS) , and  the  Surveillance  Towed 
Array  Sonar  System  (SURTASS) : -y 


Noise  Floor  Characteristics 

ASEPS  Evalvation;  Coherent  Multi  Array  Processing 
Horizontal  Directionality  of  Ambient  Noise; 

Towed  Array  Performance  , 

The  objective  of  the  Noise  Floor  Verification  Experiment  is  to 
verify  the  noise  floor  concept.  The  noise  floor  is  a depth  below 
. which  distance  shipping  noise  is  significantly  reduced  and  where 
short  range  acoustic  sensors  could  attain  a substantial  perfor- 
in mance  gain.  The  objectives  of  the  second  experiment  relates  to 

the  evaluation  of  the  propagation  model  component  of  the  Automated  \ 
Signal  Excess  Prediction  System  (ASEPS)  and  the  evaluation  of  | 

coherent  multi  array  processing  algorithms.  The  third  experi- 
ment addresses  the  measurement  of  horizontal  directionality  of 
ambient  noise  and  the  assessment  of  towed  array  performance  as 
it  relates  to  narrow  beam  noise  threshold  and  variability. 


(C)  CHURCH  OPAL  operations  were  centered  along  the  CHURCH 

ANCHOR  baseline  at  143*30 *W.  M/V  SEISMIC  EXPLORER  deployed  the 
LAMBDA  array  for  signal  propagation,  coherence,  beam  noise,  and 
noise  directionality  measurements  in  the  vicinity  of  Sites  XI, 

X2,  XA,  and  XB.  MA  AMERICAN  DELTA  II  deployed  Vibroseis  CW 
projectors  for  propagation  and  coherence  measurements  at  selected 
locations  and  along  tracks  associated  with  Sites  VI  to  Vl7a. 

R/V  MOANA  WAVE  made  a series  of  environmental  measurements  and 
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deployed  ACODACs  at  Sites  A1-A3  for  the  purpose  of  making  noise 
and  signal  propagation  measurements  at  a variety  of  depths.  The 
ACODAC  at  Site  A1  was  recovered;  the  ACODACs  at  Sites  2 and  A3 
ware  not.  In  addition  to  deploying  the  DELTA  array  for  noise 
measurements  at  Site  XB  on  one  occasion,  MOANA  WAVE  also  deployed 
the  backup  HX-231F  CW  projector  at  selected  locations  and  along 
tracks  associated  with  Sites  V13-Vl9a.  An  NRL  BP-3A  aircraft 
dropped  Mark  64  SUS  in  the  vicinity  of  the  LAMBDA  array  at  Site 
XB  on  four  towed  array  deformation  flights.  A VXN-8  RP-3A  air- 
craft made  12  environmental  (AXBT  and  ART)  measurement  and  ship- 
ping surveillance  flights.  Each  of  two  (2)  COMPATWINGSPAC  P-3C 
aircraft  made  two  (2)  shipping  surveillance  flights.  WARF  and 
SEA  ECHO  OTH  radars  conducted  surface  shipping  surveillance  opera- 
tions in  limited  areas.  Several  shore  stations  monitored  selected 
Exercise  events. 


(C)  This  progreim  is  designed  to  provide  a basis  for  manage- 

ment decisions  concerning  the  development  of  undersea  surveillance 
systems.  The  performance  of  these  systems  depends  on  the  acoustic 
environment  in  which  they  operate.  To  evaluate  the  available 
options  it  is  necessary  to  obtain  a quantitative  description  of 
these  aspects  of  the  acoustic  environment  that  affect  system  per- 
formance. The  CHURCH  OPAL  Exercise  drew  deliberately  on  the 
1973  CHURCH  ANCHOR  Exercise  data  base  and  thereby  expanded 
the  data  base  and  increased  acoustic  prediction  capability  in  the 
Northeast  Pacific  Ocean  Basin. 


(U)  A Data  Analysis  Plan* distributed  by  LRAPP  provided  speci- 

fic direction  to  participants  for  the  preliminary  phase  of  the 
analysis.  A limited  data  base  was  defined,  interpretive  analyti- 
cal techniques  were  described,  priorities  for  reduction  and  analy- 
sis of  the  data  were  enumerated,  and  schedules  and  responsibilities 
leading  to  publication  of  preliminary  analysis  results  in  March 
1976  were  established.  Budgetary  limitations  have  had  a 


LRAPP,  October  1975,  CHURCH  OPAL  Data  Analysis  Plan  (U)  (SECRET) , 
prepared  by  Xonics,  Inc. 
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significant  impact  on  this  liinited  analysis  program  and,  as  a 
consequence,  original  schedules  have  been  extensively  revised. 
Currently,  a revised  Data  Analysis  Plan  is  being  prepared  for  sub 
sequent  analysis  and  reporting  employing  a broadened  data  base.. 
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1.  (S)  INTRODUCTION  (U) 

(C)  This  report  describes  the  results  of  surveillance  of  sur>^ 

face  shipping  by  aircraft  (A/C) r comparisons  of  the  aircraft  sur- 
veillance of  shipping  with  historical  shipping  data,  and  an  in- 
vestigation of  the  feasibility  of  employing  over  the  horizon  (OTH) 
radar  for  surveillance  of  surface  shipping  during  the  CHURCH  OPAL 
exercise  in  the  Northeast  Pacific.  The  shipping  surveillance  por- 
tion of  the  exercise  was  conducted  during  September  1975,  using 
fleet  aircraft,  NAVOCEANO  VXN-8  , the  SEA  ECHO,  and  WARE  OTH 
radar  systems. 

(S)  The  objectives  of  the  CHURCH  OPAL  Exf:rcise  as  related  to 

surveillance  of  shipping  were  three-fold:  First,  to  obtain  the 

nearby  shipping  field  concurrent  with  LAMBDA  and  DELTA  horizontal 
directionality  measur^ents,  and  to  determine  the  ships  on  LAMBDA 
beams  during  beam  noise  measurements;  second,  the  comparison  of 
observed  CHURCH  OPAL  shipping  distributions  with  theoretical 
models  and  historical  shipping  fields  from  other  data  sources; 
third,  the  evaluation  of  OTH  radar  systems  to  provide  shipping  dis- 
tributions of  adequate  quality  for  use  in  theoretical  models  of 
ship  distributions  and  ambient  noise  prediction. 

(U)  The  report  is  organized  as  follows:  Sections  2 and  3 

present  the  procedure  for  reducing  and  analyzing  respectively, 
the  raw  data  obtained  from  A/C  surveillance.  The  raw  data  base 
consists  of  commentary  on  radar  and  navigation  performance,  A/C 
tracks  and  surface  ship  positions  given  by  radar  range  and  bear- 
ing from  the  A/C.  Appendix  A presents  A/C  tracks  and  ship  con- 
tacts for  the  12  flights  and  8 days  of  coverage. 

(C)  Section  4 outlines  the  procedure  used  to  generate  the 

Instantaneous  discrete  shipping  field.  A correction  procedure 
was  applied  to  some  of  the  data  which  allowed  for  a more  accurate 
reconstruction  of  ship  positions  about  the  array.  Appendix  B 
presents  the  reconstructed  shipping  fields  for  September  14  and 
16.  Individual  ship  contacts  are  given  by  range  and  bearing  from 
the  array.  Ship  positions  are  also  given  by  their  latitude  and 
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longitude:.  The  uncertainty  in  ship  positions  as  determined  by 
A/C  is  considervsd  in  Section  4. 

(0)  A comparison  between  model  and  historical  shipping  dis- 
tributions with  the  CRUltCH  OPAL  data  is  given  in  Section  5.  9%e 

reduced  CO  shipping  ddta  is  given  in  Appendix  C and  a quantita- 
tive cemparison  with  historical  data  somnarised  in  Appendix  D. 

(U)  Section  6 presents  a comparison  of  OTH  and  aircraft  ship 

distributions.  An  overall  comparison  based  upon  shipping  den- 
sities is  made  for  all  areas  for  the  entire  exercise  period,  and 
a detailed  comparison  for  September  14  between  WARE  data  and  A/C 
data  is  made  for  WARF  Area  2.  A detailed  comparative  analysici 
between  WARF  and  A/C  detections  was  possible  only  on  the  14th  of 
September  due  to  financial  constraints.  Nevertheless,  the  com- 
parison, indicated  in  Figure  €-2,  is  representative  of  the 
problems  which  are  encountered  in  attempting  to  make  a detailed 
comparison.  The  reports  of  SRI  and  NRL  on  the  operations  of 
their  WARF  and  SEA  ECHO  radars  respectively,  is  given  in  Appendix 
E. 
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2.  (C)  DATA  REDUCTION  FOR  SHIPPING  DENSITIES  (U) 

(U)  The  location  of  observed  surface  vessels  and  the  tracks 

flown  by  the  surveillance  aircraft  plus  estimates  of  the  per for ~ 
mance  of  the  A/C  radar  system  at  each  point  where  ships  were  ob" 
served  by  radar  along  the  flight  track  represents  the  raw  data 
base.  From  these  raw  data  three  numbers  were  derived  for  each 
five  by  five  degree  square  within  A/C  radar  range  for  each  leg  of 
each  flight.  The  first  is  the  number  of  contacts  observed  in  a 
given  square  within  radar  range  for  that  leg»  the  second  is  the 
portion  of  the  given  square  within  the  estimated  radar  range  of 
the  flight  path,  end  the  third  is  ship  density  for  the  given  square. 
A ship  density  is  not  calculated  if  leas  than  half  the  square  was 
surveyed. 

(U)  Overlapping  areas  created  at  A/C  track  turn  points  were 

apportioned  approximately  equally  between  the  two  legs  in  order 
not  to  be  counted  twice.  From  observeri^  comments  concerning  the 
generally  good  radar  operation,  a standard  radar  coverage  radius 
of  75  nm  was  established  and  was  used  except  in  those  cases  where 
it  was  obvious  from  the  observer,  navigation  or  contact  logs  that 
a reduced  or  greater  radar  range  was  in  effect.  As  a check  on 
radar  operations,  a straightforward  statistical  test  may  be  em- 
ployed. The  aircraft  track  is  generally  a sequence  of  rhumb  line 
segments  which  crisscross  the  area  of  Interest.  Although  a leg 
may  parallel  a shipping  route  for  a short  distance,  when  one  con- 
siders a statistic  involving  the  entire  track,  these  local  fluc- 
tuations are  smoothed  out  by  the  law  of  large  numbers.  For  each 
of  the  N ship  contact  observed  on  a flight,  there  is  a 

distance  r^,  the  range  of  the  ship  from  the  flight  track.  If 
the  radar  has  good  coverage  to  75  nm  (or  more) , then  the  set  of 
all  ranges  observed  on  the  flight  which  were  less  than  75  nm 
should  be  approximately  uniformly  distributed,  l.e.,  the  set 

{rj^  s 1 <.  i <.  N,  ^ 

should  be  representative  of  a draw  from  a*a  (0,  75]  distribution. 
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If  the  distribution  is  not  uniform  but  decreases  with  range, 
then  it  can  be  assumed  that  surface  vessels  within  the  coverage 
radius  are  being  missed  and  a resultant  bias  is  being  intro- 
duced. A determination  that  these  distributions  were  approxi- 
mately uniform  preceded  the  estaJalishment  of  the  expected  ra- 
dar coverage  range. 

(U)  For  each  exercise  day  the  ntodser  of  contacts  and  the 

percent  coverage  for  each  five  by  five  degree  square  were  summed 
over  all  legs  of  all  flights,  e.g.,  if  a flight  covered  75%  of 
a specific  square  on  both  an  outbound  and  inbound  leg  then  per- 
centage coverage  due  to  that  flight  would  be  150%  and  the  num- 
ber of  contacts  would  be  the  sum  of  the  two  observations.  It 
is  recognized  that  any  multiple  coverage  thus  created  does  not 
necessarily  provide  completely  independent  observations  since 
the  relaxation  time  is  a function  of  the  average  ship  transit 
time  through  the  squa.:e;  however,  this  procedure  does  not  intro- 
duce any  bias  and  the  variance  of  the  density  estimator  is  re- 
duced by  any  deviation  from  a correlation  of  1.0  that  might 
exist  between  observations. 

(C)  Ship  positions  and  their  range  and  bearing  from  the 

LAMBDA  array  on  the  dates  of  14  and  16  September  are  given  in 
Appendix  B.  Array  location  and  tiroes  are  given  in  Table  Bl.  Due 
to  A/C  navigation  limitations,  it  was  decided  that  array  location 
be  determined  from  the  M/V  SEISMIC  EXPLORER  Satellite  Navigation 
Log.  The  actual  ship  positions  to  be  used  in  determining  LAMBDA 
beam  noise  threshold  levels  are  given  in  Tables  B2  and  B3.  These 
ship  positions  were  checked  against  uhe  Naval  Oceanographic  Office 
(NOO)  exercise  reconstruction,  and  it  was  found  that  the  positions 
were  in  agreement. 

(U)  Ship  range  was  computed  for  each  surface  contact  using 

equation  (1)  for  the  great  circle  distance  between  the  array  and 
each  ship  location.  Bearing  as  measured  from  the  array  was  com- 
puted using  equation  (2),  where  latitude  and  longitude  of  the 
array  is  given  by  (0  , ) and  for  the  surface  ship  (C  , ^ ). 


D(nm)  » 60  cos'^Esin  0^  sin  6^  + cos  6^  cos  0^  cos(^^  - ^1^)3  d) 

sin  0g  - sin  0^  cos (D/60) 

® ■ • ET5-|d/80)  Sos-B"  ' 

Cl 


B ; sin(4»g  - < 0 

360  - B 1 sin((J)  - ♦ ) > 0 

9 A 


(2) 


(U)  The  sum  of  the  contacts  for  each  five  by  five  degree 

square  is  divided  by  the  sum  of  coverage  over  that  square  to  give 
the  estimated  ship  density  for  the  given  exercise  day.  The  re- 
sults are  contained  in  Appendix  C which  presents,  for  each 

measuronent  day: 


(1)  the  sum  of  the  contacts 

(2)  the  sum  of  the  coverage 

(3)  the  estimated  shipping  density  (i.e., 
the  quotient  (1)/(2)J  for  all  x 5® 
squares  having  coverage  of  50%  and 
greater . 


We  note  that  the  estimates  are  given  for  shipping  densities  — not 
the  ntimber  of  ships  in  a square  which  are  customarily  estimated  in 
historical  shipping  "density"  charts.  The  two  terms  are  identi- 
cal except  when  land  area  is  Included  within  the  square  and  for 
such  a square  an  adjustment  must  be  made  before  a comparison  with 
historical  data  can  be  undertaken. 

(C)  The  actual  noise  directionality  is  computed  from  a di- 

rectional ambient  noise  model.  These  models,  of  which  there  are 
a small  number,  require  as  Inputs  shipping  densities  in  some  cases, 
and  actual  ship  positions  in  other  cases.  In  those  models  which 
require  shipping  densities,  the  data  is  provided  in  1"  x 1"  squares. 
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That  re<sulrea  that  the  5**  x 5*  square  data  be  broken  down  to 
1*  X 1*  square  data.  This  Is  done  using  an  exponential  weighted 
smoothing  algorithm.*  The  models  then  accept  this  1**  x 1*  square 
data  derived  from  the  baaic  5*  x 5*  square  data  and  breaks  it 
down  into  even  mailer  areas,  vis.,  10*  x 10*  square.  This  is 
obtained  from  the  1*  x 1*  square  data  generally  by  uniformly  dis" 
tributing  the  densities  into  the  36  - 10*  x 10*  squares  in  each 
1”  X 1*  square.  Ccx^iarisons  in  determining  the  validity  of  ship- 
ping surveillance  data  is  almost  always  made  by  comparing  5**  x 5^ 
square,  or  even  larger  areas.  The  smaller  number  of  ships  in  any 
1^x1*  square  generally  is  too  variable  to  allow  reasonable 
comparisons  with  various  data  and  these  larger  areas  are  almost 
always  employed.  There  are  exceptions  to  this  procedure,  but  only 
in  very  special  circumstances  such  as  in  narrow  seas,  straits,  or 
similar  confined  ocean  areas. 

(C)  The  use  of  shipping  densities  as  opposed  to  real  or  spe- 

cific ship  positions  in  the  ambient  noise  xnodels  i%  to  a great 
extent,  dependent  upon  the  system  which  is  being  studied . If  a 
system  with  broad  beams  is  being  considered,  then  there  is  always 
a large  number  of  ships  seen  within  that  beeuh.  Therefore,  if  a 
few  ships  are  excluded  from  consideration,  the  effect  on  the  beam 
noise  level  is  quite  small.  At  the  low  frequencies  which  are  being 
considered,  the  propagation  loss  is  not  of  great  importance  to 
the  ambient  noise  field,  and  it  may  be  fairly  argued  that  the  inte- 
gration of  the  shipping  density  over  the  area  covered  by  the  beam 
is  the  major  effect  on  the  beam  noise  level.  Comparison  of  mea- 
sured with  predicted  beam  noise  levels  Indicate  the  soundness  and 

*(L.P.  Solomon,  "CHURCH  ANCHOR:  AIRCRAFT  SURVEILLANCE  OF  SHIPPING" 
PSI  TR-004002,  Pages  13-18). 
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validity  of  this  approach.  A direct  result  of  this  observation  is 
the  recognition  that  the  accuracy  constraints  on  the  shipping  den- 
sity is  5*  X 5**  square  is  much  less  stringent  than  on  a more  com- 
pact density  specification.  Our  experience  with  the  5**  x 5^  his- 
torical shipping  fields  as  modified  by  at-sea  exercise  data  indi- 
cates its  dependability  and  utility  when  used  by  directional  ambi- 
ent noise  models  in  considering  broad  beam  systems. 


UnCLASS3B1£D 


3.  (U)  STATISTICS  OF  THE  SHIPPING  DENSITIES  (0) 


(U) 


For  each  5**  x 5*  square  let  be  the  number  of  ships 

JL 


.th 


observed  on  the  day  (i  « 1,2,...K)  and  K the  number  of  days 

of  coverage  for  a given  square • Define  to  be  the  proportion 

of  the  area  covered  on  the  i^  day.  Tf  p.  0 then  an  estimate 

i th 

of  the  density  of  ships  in  the  area  on  the  i—  day  is  given  by 


\ - r.^/p^ 


(3) 


The  daily  weighted  ship  densities  are  given  in  Appendix  C,  Figures 
Cl  through  C8.  For  later  convenience,  define: 

K 

» = Pi  (4a) 

i=l 
K 


(4b) 


(U) 


i=l 

The  weighted  average  has  the  form 

K 


-i-  2 V 

m 1 


(5) 


i«l 


41 

1 


where  w^  and  m are  constants  (the  constant  m is  introduced  so  that 
the  Wj^  may  have  simpler  form) . In  order  to  attach  the  same  impor-  | 
tance  to  each  square  mile  which  was  covered,  the  weight  must  be  | 

proportional  to  the  amount  of  area  covered , i . e . , w^^  - p^ . The  J 

constant  m is  chosen  such  that  d is  an  unbiased  estimator.  That  1 

is,  if  d^  are  independent  samples  of  a random  variable  (called  1 

"density")  having  mean  y,  then  m is  chosen  such  that  the  expected  I 
value  of  d is  y,  i.e.,  I 


E[d]  * y (6) 

which  implies  m * s.  But  w^d^  * ^i^i  “ *^i 

K 


(U)  For  similar  reasons,  a vreighted  estimator  of  variance  is 

desirable.  This  estimator  has  the  form 


S Wj.  (d^  “ d)^  (8) 

i-1  ^ 

Letting  gives  a system  consistent  with  equation  (7).  This 

system  favors  the  larger  area  coverage,  but  does  not  place  inordi- 
nate emphasis  on  the  days  with  maximum  coverage.  Choosing  ro  such 
that  9^  is  an  unbiased  estimator  gives 


e2  = 


s - s 


< E n^/Pi  “ 


2 i»l 


where  n^/p^  = 0 when  n^  = p^^  = 0. 

(U)  Ship  density  statistics  for  5**  x S'*  squares  are  provided 

if  there  was  50%  coverage  or  more  for  the  given  square.  The  num- 
ber of  days  K,  for  which  A/C  coverage  met  the  above  criterion  for 
the  5®  X 5*  squares  is  given  in  Appendix  D,  Figure  Dl.  We  note 
that  there  is  a variation  in  K between  squares,  which  depends  on 
the  number  of  flights,  radar  range,  and  the  extent  of  overlapping 
coverage  in  the  given  area.  Weighted  average  densities  and 
standard  deviations  for  these  squares  are  given  in  Figure  D2  and 
Figure  D3.  Since  the  sums  in  Equations  (7)  and  (9)  are  over  K,  our 
confidence  in  these  statistics  must  vary  accordingly. 

It  is  of  interest  to  compare  the  CHURCH  OPAL  (CO)  averaged 
distributions  with  predicted  model  results  and  other  historical 
files.  Such  a comparison  may  be  used  to  evaluate  the  model  pre- 
dictions, note  any  temporal  variability  in  the  shipping  fields  and 
to  up-date  the  model  shipping  routes  if  necessary.  The  averaged 
shipping  densities  obtained  in  September  of  1973  during  the  CHURCH 
ANCHOR  (CA)  exercise  together  with  the  RMS  and  Automated  Marine 


International  (AHI)  model  distribution  are  given  for  5*  x 5* 
squares  in  Figure  D3.  The  ratio  of  CA,  RMS  and  AHI  densities  to 
CO  is  given  in  Figure  D4.  Analysis  of  these  data  is  given  in 
Section  5. 
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4.  (S)  AIRCRAFT  PERFORMANCE  AMD  LIMITATIONS  (U) 

(C)  Radar  coverage  of  surface  vessels  was  provided  on  eleven 

days  during  the  CHURCH  OPAL  Exercise.  On  each  of  two  high  priori- 
ty LAMBDA  beam  noise  measurement  days  (14  and  16  September) , three 
flights  covered  the  surrounding  area.  OTH  radar  provided  overlap- 
ping coverage  on  the  14th  for  the  region  labeled  Area  2 in  Figure  El. 

(C)  Appendix  A gives  the  actual  tracks  for  each  flight  and 

Table  A1  gives  the  schedule  of  flights  for  each  exercise  day.  Re- 
liable data  were  not  obtained  on  Flight  PB  on  the  14th  due  to  radar 
malfunction.  The  flight  was  Intended  to  provide  far-field  shipping 
data  for  the  region  18-22 ®N,  137-158 ®W.  Fortunately,  excellent  ra- 
dar conditions  prevailed,  and  flight  PA  was  able  to  provide  partial 
coverage  in  this  region.  Ship  contacts  were  obtained  in  this  re- 
gion which  are  over  950  nm  from  the  array.  Flights  on  September  26, 
29  and  October  1 were  at  very  high  and  low  altitudes,  and  not  pri- 
marily intended  for  shipping  surveillance/  but  for  airborne  radia- 
tion thermometry  measurements.  Consequently,  no  statistics  are 
provided  on  these  dates;  the  other  eight  days  were  extensively 
analysed . 

<U)  Flight  Commanders  were  requested  to  fly  at  altitudes  con- 

sistent with  both  good  radar  coverage  and  low  clutter  rate  for  the 
prevailing  conditions.  A scientific  observer  on  each  flight  was 
responsible  for  the  collection  of  the  radar  contact  logs  and  navi- 
gation logs  as  well  as  for  a personal  log  in  which  he  was  to  docu- 
ment any  deviations  from  the  scheduled  operations,  e.g.  radar  mal- 
functions, changes  in  flight  plan,  bad  weather  conditions,  etc. 

(C)  For  each  flight  an  overall  mission  effectiveness  is  given 

in  Table  Al.  These  evaluations  are  assessments  over  the  entire 
duration  of  the  flight.  Weather  conditions  may  have  varied  dra- 
matically over  certain  portions  of  the  track.  Therefore,  these 
performance  estimates  provide  a measure  of  the  mean  effectiveness 
for  a given  flight.  Accordingly,  we  judge  overall  A/C  surveillance 
for  these  eight  days  of  coverage  to  have  been  very  good.  The  crews 
of  both  VXN“8  and  VP-1  were  excellent  and  their  performance  is  re- 
sponsible for  the  overall  success  of  the  missions.  Of  particular 
merit  was  the  outstanding  performance  of  personnel  at  the  Tactical 


Support  Center,  Barbers  Point. 


(S)  In  the  following  we  delineate  the  procedure  used  to  gen- 

erate the  Array  and  Ship  Position  Tables  Bl,  B2  and  B3.  The  ref- 
erence times  given  in  these  tables  for  14  and  16  September  are  the 
times  A/C  contact  was  made  of  the  M/V  SEISMIC  EXPLORER  fSE) . Using 
data  from  the  SE  Satellite  Navigation  Log  (SATNAVLOG)  for  ship 
position  and  speed,  the  position  of  SE  was  determined  at  the  ref- 
erence times  by  dead  reckoning  (DR) . Ship  logs  showed  that  SE  did 
not  change  course  or  speed  during  the  time  interval  over  which 
the  DR  computation  was  made  so  that  we  have  high  confidence  in  the 
SE  positions  given  in  Table  Bl  which  for  2336A  14  September  was 


31.6378®  N,  143.0780®  W. 

The  position  of  SE  determined  by  VXN-8  flight  TE  on  the  14th  from 
a range  of  75  nm  using  APS-80  radar  at  the  same  time  was 

31.6055®  N,  142.9521®  W 


Thus  the  two  positions  differ  in  latitude  by  AL  = 0.0323®  and 
longitude  by  6X  0.1259®.  The  difference  in  range  6r  between 
the  two  positions  {using  Equation  (1)]  is  6.72  nm.  Assuming  no 
uncertainty  in  dead  reckoning  the  SE  position  from  Satellite  Navi- 
gation data,  this  distance  can  be  thought  of  primarily  as  the  error 
in  the  aircraft's  track.  The  array,  of  course,  is  used  in  the 
post-exercise  reconstruction  as  a tie  point. 

(S)  On  14  September  VXN-8  navigation  was  dependent  upon  a 

ccsnblnation  of  LORAN-C  (which  was  noted  to  be  very  poor  for  the 
above  contact),  celestial  and  radio  fixes.  On  16  September  VXN-8 
employed  the  LITTON-51  inertial  navigation  system.  Contact  of 
the  SE  was  logged  at  a range  of  19  nm.  An  analogous  computation 
gave  differences  in  latitude  6L  = 0.0591®  and  longitude 
6X  “ -0.0010®.  This  positional  discrepancy  represents  a range 
difference  of  6r  = 3.54  nm  between  the  two  locations. 

(C)  For  all  surface  contacts  logged  by  VXN-8  on  the  above 

dates  the  corresponding  latitude  and  longitude  differences,  5l, 
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6X  were  added  to  all  A/C  ship  positions.  This  correction  is 
only  valid  for  VXN**8  contacts;  ship  positions  reported  by  fleet 
A/C  were  not  altered.  The  shift  in  position  of  the  VXN-8  con” 
tacts  allowed  a more  accurate  reconstruction  of  the  shipping. 

(U)  Tables  b2  and  B3  are  intended  to  provide  the  Instantaneous 

discrete  shipping  fields  for  times  2336Z  14  September  and  2259Z 
16  September  1975 , respectively.  The  positions  of  all  surface  con- 
tacts were  not  determined  simultane<^usly , but  over  a period  of 
several  hours  during  which  the  shipping  field  did  not  remain  sta- 
tionary. Since  surface  ship  speeds  and  bearings  were  not  obtained, 
precise  reconstruction  of  the  discrete  shipping  field  is  not  possible. 
Estimates  however,  can  be  made  that  provide  an  upper  bound  to  the 
uncertainty  introduced  into  the  range  and  bearing  of  ship  positions 
as  measured  from  the  array. 

(U)  Uncertainty  in  position  of  a surface  contact  with  respect 

to  the  array  will  be  introduced  into  Table  B2  and  B3  if  there 
occurred  relative  motion  between  the  contact  and  array  over  the 
time  interval  between  A/C  recording  of  the  contact  and  the  reference 
time  given  in  the  tables.  If  the  relative  motion  was  transverse, 
we  approximate  the  upper  limit  in  bearing  uncertainty  6b  by 

6B  = v(At)/x  (10) 

where  v is  the  surface  ship  speed.  At  the  time  interval  between 
A/C  determination  of  contact  and  the  reference  time,  and  x the 
radius  vector  from  array  to  surface  contact.  Had  the  relative  mo- 
tion been  entirely  radial,  the  range  uncertainty  6R  is 

6R  = v(At)  (11) 

(C)  As  representative  exeunples  we  take  from  Table  B2  entries 

1,  16  and  20  and  compute  6b,  6r  for  these  cases.  The  results  are 
summarized  in  Table  4.1  Entrees  1,  20  and  16  are  for  ship  contacts 
in  the  far,  middle  and  near  shipping  fields  respectively.  They 
are  919.1,  130.6  and  14.7  nm  from  the  array.  Assuming  these  con- 
tacts were  moving  transverse  to  the  SE  and  therefore  parallel  to 
the  array,  the  maximum  bearing  uncertainties  are  3.1,  3.5  and  0.78" 
respectively.  Thus  we  could  write  for  entry  16  a bearing  of 
38.0769  i 0.7802"  (bearing  has  been  measured  from  0"  T. 
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Had  the  motion  been  entirely  longitudinal , the  range  would  be  given 
by  14.6680  it  0.2000  nm.  These  examples  illustrate  the  gross  bear- 
ing and  range  resolution  uncertalnticus  that  may  be  expected  and 
the  possible  limitations  of  the  discrete  shipping  fields  given  in 
Tables  B2  and  B3. 
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Table  4.1  (C) 

Maximum  Range  and  Bearing  Uncertalntlea  due  to  Ship  Movement 
for  Table  B2  Entries  1,  16  artd  20*  (tJ) 

(Reference  time  142336  Z) 


Table  B.  2 
Entry 

Range  from 
Array  (nm) 

Time  Logged 
(DTG) 

6R'*' 

(nm) 

6B++ 

(deg.) 

1 

919.0800 

141931 

48.96 

3.0539 

16 

14.6880 

142337 

0.200 

0.7802 

20 

130.6212 

150016 



8.000 

3.5093 

* Assiunes  surface  ship  speed  of  12  knots  | 

”*■  Assumes  relative  motion  is  entirely  radial  i 

Assumes  relative  motion  is  entirely  transverse  | 
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5.  (O)  COMPARISON  WITH  MODEL  AHD  HISTORICAI.  i3H'»J>FIWG  DBMSITiaS  (O) 

(U)  In  this  section  we  examine  the  variability  la  shi]^.>ir9 

densities  for  the  CHURCH  OPAL  Exercise  are  examined.  TL.>  «A''..ghted 
CO  densities  are  compared  with  historical  CHURCH  ANCHOR,  *'.adel  RMS 
and  AMI  shipping  fields. 

(U)  ' Examination  of  Figures  C1-C8  shows  that  the  individual 
CO  weighted  ship  densities  [d^^  of  Equation  (3)]  within  a given 
5*  X 5**  square  may  be  highly  variable.  Furthermore,  it  seems  that 
this  variability  can  be  associated  wirh  coastal  and  open  ocean  re- 
gions alike.  Prom  Figures  D2  and  D3  for  the  weighted  average  den- 
sity and  standard  deviation  (SD) , we  note  a density  variation  in 
range  of  over  100%  from  squares  of  low  and  high  standard  deviation.* 
Therefore,  it  would  appear  that  in  the  near-field,  i.e.,  within 
a 5*  X 5**  square,  the  shipping  distribution  may  be  well  respre- 
sented  as  a stochastic  process. 

(U)  The  overall  ship  densities  computed  during  CHURCH  ANCHOR 

and  CHURCH  OPAL  show  remarkable  agreement.  To  study  the  results 
in  detail  we  consider  two  regions,  separated  by  the  35*N  latitude. 
Above  this  line  the  squares  with  common  (i.e.,  CA  and  CO)  coverage 
generally  had  greater  coverage  during  the  CHURCH  ANCHOR  exercise 
than  during  the  CHURCH  OPAL.  Conversely,  those  south  of  35®N 
had  far  better  coverage  during  CHURCH  OPAL  (indeed,  many  of  the 
lower  squares  shown  in  Figure  D4  had  no  coverage  during  CHURCH 
ANCHOR) . In  the  southern  regioif  the  CHURCH  OPAL  estimates  gave  a 
total  of  57  ships,  as  compared  with  62  for  CHURCH  ANCHOR.  However, 

* Write  the  average  weighted  density  as  d 1 SD  and  define  per  cent 
variation  for  a square  as  (SD/d)  x 100.  Using  Figures  D2  and  D3 
for  squares  20-25®N,  140-145®W;  30-35*N,  145-150®W  the  per  cent 
variation  in  density  is  5.69  and  106.98  respectively.  Therefore 
the  range  of  the  variation  (in  percent)  between  these  squares  is 
101.29. 

“^To  compare  CO  and  CA  densities  throughout  the  CO  region,  it  was 
necessary  to  use  CA  densities  extended  to  include  roughly  the 
area  20-30®N,  140-160®W.  The  algorithm  used  to  extend  these 
densities  to  areas  not  covered  during  the  CA  exercise  is  given  by 
L.P.  Solomon  (PSl  Report  TR-004002) . Figure  C3  and  Dl  of  this 
reference  were  used  to  give  the  AMI  and  CA  extended  densities. 
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th«re  was  a significant  difference  in  the  distrlbation  of  these 
ships  between  the  squares.  The  shipping  in  the  neighborhood  of 
Hawaii  and  Northwest  of  the  islands  Is  distributed  over  a sonie''  ^ 
what  broader  region  than  either  the  CHURCH  ANCHOR  or  WIS  fields 
show  (both  of  these  distributions  show  rather  narrow  shipping  Iwes 
from  Hawaii  to  the  mainland) . It  must  be  remembered  that  the 
shape  of  the  CHURCH  ANCHOR  distribution  (but  not  the  total  numbers) 
in  this  area  was  determined  from  the  AMI  data,  since  there  was  no 
shipping  surveillance  coverage  during  that  exercise  in  the  immedi- 
ate vicinity  of  Hawaii.  This  data  tends  to  confirm  the  thesis 
that  shipping  lanes  in  the  open  ocean  are  not  as  “narrow"  as  is 
generally  believed*  Other  exercises  have  also  Illustrated  this 
phenomenon.  In  particular  CHURCH  ANCHOR  showed  that  the  West  Coast- 
Japan  route  was  much  broader  and  dyneunic  than  Indicated  by  the 
synthesized  shipping  fields. 

(U)  In  the  region  north  of  35"N,  CHURCH  OPAL  gives  a total 

of  56.95  ships  versus  56.12  for  CHURCH  ANCHOR.  When  comparing 
densities  (Figure  D4)  it  is  seen  that  in  the  three  coastal  squares 
the  CHURCH  OPAL  figures  are  lower  than  those  of  CHURCH  ANCHOR. 

These  in  turn  are  lower  than  the  RMS  numbers  for  these  squares, 
vhich  are  (generally)  less  than  the  AMI  numbers.  The  CHURCH  ANCHOR 
data  base  was  extensive  in  these  coastal  squares,  but  CHURCH  OPAL 
was  quite  light;  thus  the  small  sample  does  not  justify  a conclu- 
sion that  coastal  traffic  is  changing. 

(U)  In  conclusion,  the  shipping  fields  from  the  CHURCH  OPAL 

exercise  (Sept.  1975)  agree  remarkably  well  with  those  of  CHURCH 
ANCHOR  (Sept.  1973).  The  CHURCH  OPAL  data  in  the  vicinity  of 
Hawaii  indicates  the  actual  nature  of  the  Hawaii-roainland  shipping 
lanes,  superseding  the  previous  synthesized  distributions. 
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4,  (S)  COe^ARISDN  (Hr  OTH  AMD  AXRCltAFT  SHIP  DISTRIBUTIONS  <U) 

(C)  The  location  of  the  CHURCH  OPAL  Exercise  Area  allowed 

for  surveillance  of  surface  shipping  by  the  SEA  ECHO  dver^Tlie' 
Horizon  (OTH)  radar  system  at  San  Clemente  island  and  the  Wide 
Aperture  Research  Facility  (WARF)  at  Los  Banos » California.  Sur- 
veillance of  the  relatively  large  ocean  areas  (5^  x 5**  square)  rer 
qulred  strategic  scanning  and  subsequent  doj^ler  processing  to 
discriminate  individual  ships  from  the  sea  return.  The  detection 
criterion  is  based  on  measurement  of  the  radial  component  of  ship 
velocity,  foz^  which  below  a threshold  value,  the  ship  will  not 
be  observed.  For  the  SEA  ECHO  radar,  ships  having  a radial  veloci- 
ty of  less  than  approximately  13  knots  would  be  below  the  detec- 
tion threshold.  Thus  OTH  radar  may  underestimate  the  ship  densi- 
ty according  to  the  directional  distribution  of  the  shipping  den- 
sity. 

(C)  On  September  14,  22  and  24  there  was  mutual  coverage  of 

shipping  by  OTH  radar  and  A/C  in  the  areas  shown  in  Figure  El. 

The  WARF  OTH-B  system  provided  surveillance  on  the  14th  and  SEA 
ECHO  radar  system  on  the  22nd  and  24th.  Table  El  summarizes  the 
results  of  OTH  and  A/C  shipping  surveillance  operations  for  the 
three  areas  on  the  above  dates. 

(S)  Analysis  of  WARF+' shipping  distributions  (J.  R.  Barnum, 

Ship  Density  Determination  with  High  Resolution  Skywave  Radar 
Surveillance  (SECRET) , October  1975)  indicates  the  presence  of 
12  ships  within  WARF  Area  2.  Aircraft  coverage  in  Area  2 for 
time  comparable  with  OTH  gave  a total  of  11  surface  contacts. 

Since  A/C  radar  conditions  were  excellent  on  14  September,  the 
confidence  level  for  these  contacts  is  >90%.  Most  of  these  con- 
tacts appiwar  to  have  been  observed  on  OTH;  others  were  not. 

There  are.  10  probable  correlations.  Reliability  of  OTH  data  will 
suffer  unless  sufficient  time  for  on-line  analysis  and  verifi- 
cation of  each  possible  contact  is  provided.  A plot  of  contacts 
for  14  September  obtained  from  multiple  sources  (VP,  VXN-8,  and 


'^See  Appendix  E. 
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OTH)  is  given  in  Figure  E2.  Correlation  of  multiple  source  con- 
tacts is  a donanding  exercise  in  judgment. 


(C)  S£A  ECHO  coverage*'  (R.  W.  Bogle,  OTH  Radar  Ship  Distri- 

bution Measurements;  September  1975)  on  the  22nd  and  24th  is  shown 
in  Figures  E3  and  £4.  The  three  OTH  contacts  ware  dead  reckoned 
3-4  hours  assuming  12  knot  ship  speed  and  corresponded  to  3 of 
VXN-8's  contacts.  The  two  outside  0TB  contacts  were  dead  reckoned 
to  points  outside  of  VXN-8  radar  range.  The  remaining  four  VXN-8 
contacts  were  not  detected  by  OTH,  and  they  were  too  far  frcxn  the 
edge  of  the  area  to  have  moved  "into  the  region"  during  the  time 
difference.  Weather  for  VXN-8  was  quite  poor. 

(C)  On  September  24  OTH  found  4 contacts  in  the  area  — 4 others 

immediately  outside.  VXN-8  found  3 contacts  in  the  area.  Since 
the  time  difference  is  so  great  and  only  radial  velocity  was  de- 
tected, it  is  not  possible  to  make  an  accurate  comparison.  How- 
ever of  the  8 OTH  contacts,  it  could  be  that  5 were  spotted  by 
VXN-8  (3  in  the  square  and  2 outside) , and  1 moved  out  of  the  area, 
leaving  2 OTH  contacts  which  could  not  have  moved  out  of  VXN-8 's 
range  but  were  not  detected  by  them. 


Date 

Coverage  Times  (Z) 
(QTH/Aircraft) 

14  Sept. 

1639-0102 

2130-0033 

22  Sept. 

1800-1830 

2157-2312 

24  Sept. 

1430-1500 

2319-2335 

24 
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7.  (S)  CONCIiUSIONS  (U) 

(S)  The  objectives  of  the  CHURCH  OPAL  Exercise  as  related  to 

surveillance  of  shipping  were  three-fold:  first,  to  obtain  the 

nearby  shipping  field  concurrent  with  LAMBDA  and  DELTA  horizontal 
directionality  measurements,  and  to  determine  the  ships  on  LAMBDA 
beams  during  beam  noise  measurements;  second,  the  ccnnparison  of 
observed  CHURCH  OPAL  shipping  distributions  with  theoretical 
models  and  historical  shipping  fields  front  other  data  sources; 
third,  the  evaluation  of  OTH  radar  systems  to  provide  shipping  dis- 
tributions of  adequate  quality  for  use  in  theoretical  models  of 
ship  distributions  and  ambient  noise  prediction. 

(S)  The  conclusions  which  can  be  drawn  from  this  study  can 

be  divided  into  three  major  areas:  data  obtained  from  aircraft 

determination  of  the  discrete  shipping  fields;  data  generated  from 
the  statistics  of  the  shipping  distributions  as  contrasted  with 
observed,  historical  and  model  predictions;  and  data  used  for  com- 
parative assessment  of  OTH  radar. 

(1)  On  the  two  high  priority  beam  noise  measurement 
days,  errors  in  ship  positions  were  small  with 
known  bounds. 

(2)  Operational  conditions  were  such  that  overall 
surveillance  on  the  eight  days  of  aircraft 
coverage  was  good. 

(3)  The  ship  densities  from  the  CHURCH  OPAL 
Exercise  (Sept.  1975)  are  consistent  with 
those  of  CHURCH  ANCHOR  (Sept.  1973).  Both 
CHURCH  OPAL  and  CHURCH  ANCHOR  shipping  densities 
are  in  good  agreement  with  the  RMS  distri- 
bution, except  near  the  coasts  where  greater 
variation  was  expected. 

(4)  During  the  summer,  the  historical  shipping 
density  (RMS)  in  the  Northeast  Pacific  is 
judged  to  be  sufficiently  precise  to  serve 


ev 


as  an  Input  to  present  open  ocean  low 
frequency  broad  beam  directional  axnblent 
noise  models. 

(5)  The  OTH  radar  is  a potentially  valuable 
tool  for  shipping  density  estimation, 
but  should  presently  be  utilised  only 
in  conjunction  with  aircraft  surveillance. 
There  is  a problem  of  correlation  between 
aircraft  and  OTH  contacts.  Resolution 
problems  remain  with  OTH;  however,  these 
may  be  sharply  alleviated  by  increased 
dwell  time. 
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TABLE  A1  (C)  FLIGHT  SUMMARY  (U 


Date 

Plight 

Track 

Originations 

Destination 

Mission 

Status* 

10  Sept. 

TC 

Moffett 

■ 

Barbers  Pt. 

2 

12  Sept. 

TD 

Barbers  Pt. 

Barbers  Pt. 

2 

14  Sept. 

TE 

Barbers  Pt. 

Barbers  Pt. 

1 

14  Sept. 

PA 

Barbers  Pt. 

Barbers  Pt. 

1 

14  Sept. 

PB 

Barbers  Pt. 

Barbers  Pt. 

NC 

16  Sept. 

TE 

Barbers  Pt. 

Barbers  Pt. 

1 

16  Sept. 

PA 

Barbers  Pt. 

Barbers  Pt. 

2 

16  Sept. 

PB 

Barbers  Pt. 

Barbers  Pt. 

2 

18  Sept. 

TG 

Barbers  Pt. 

Moffett 

1 

20  Sept. 

1 

TA 

Moffett 

Elmendorf 

PCI 

22  Sept. 

1 TH 

Elmendorf 

Moffett 

PC2 

24  Sept. 

1 TA 

Moffett 

Elmendorf 

1 

*General  Status  codes  for  overall  flight  effectiveness  are: 

1 = Overall  excellent  radar  and  navigation;  coverage  along  entire 

track 

2 = Overall  good  radar  and  navigation;  coverage  along  entire  track 

3 =>  Overall  poor  radar  and  navigation;  coverage  along  entire  track 

PCI  =»  Excellent  radar  coverage  along  95%  of  track;  no  coverage  along 
the  track  through  35-36®N,  135-138®W 

PC2  - Excellent  radar  coverage  along  60%  of  track;  no  coverage  over 
remainder 

NC  = No  coverage 
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FlGmtB  AIS  (C)  TRACK  iPA  14  SEPTEM8BK  W) 

iBL  1^  im  iHE.  iHft  las.  laoL  ^ tm 


DATE 


TIME 

(Zulu) 


LOCATION 

(Latitude,  Longitude) 


14  September 

2336 

31.6378*N,  143.0780»W 

16  September 

2259 

31.5452»N,  144.2988*W 

TABLE  B1  (C)  Location  of  LAMBDA  Array 

Used  to  Determine  Ship  Ranges 
and  Bearings  (C) 


. J 


i' 


latitude  ■ 22.4146 
longitude  “156. 840S 
rarigf  •'919.0800 
baariJQtg  ^36.3570 

21.8804 

156.27913 

916.7118 

233.5348 

22.4613 
156.2211 
891.1066 
235.0659 

23, 4135 
155.2076 
811.4156 
235.5446 

23.8456 
154.2561 
754.7797 
234.4990 

23.7131 
150.5433 
618.8166 
221.6388 

24.8522 
153. 0656 
666. 0836 
234.8179 

24.8425 
149.8821 
543.3826 
223.0799 

24.2905 
147.4502 
497.8763 
208.7813 

26.6072 
144.6695 
313.  1344 
195.3425 


28.0641 

146.9560 

294.4000 

224.2457. 

28.2842 

145.2608 

230.9841 

209.9708 

27.7647 

144.4480 

243.0990 

197.4218 

26.7249 

144.4448 

303.3305 

193.9904 

28.0590 

146.3160 

272.9149 

218.9409 

31.8304 

142.9004 

14.6880 

38.0769 

30.5210 

142.2978 

78.0847 

148.9036 

30.7932 
142. 1440 
69.7507 
136.3527 

30.9156 

1-43.4926 

43,2616 

206.2395 

33.8137 
143. 1620 
130.6212 
358. 1648 


33.3704 

144.7761 

134.8675 

320.8790 

33.8378 
143. 1749 
132.0920 
357.9061 

36. 1822 
143.4636 
273.3364 
356.0793 

37.6280 

143.0780 

359.4104 

0.0000 

37.6513 

144.3681 

366.3736 

350.3533 

36.8357 

145.3183 

331.0484 

341.0106 

38.5929 

144.6903 

424.7223 

349.7207 

36.9144 

145.8155 

344.4132 

337.5561 

36.5892 

147.3855 

366.0192 

325.4248 

34.5371 

146.2541 

236,0912 

318.3105 


32.7766 
149.4486 
338. 4972 
288r6165 

"33.8395 

153.0624 

520.6806 

287.3590 

25.3460 

151.4575 

580.7850 

231.5608 

23.9630 

154.2361 

749.3706 

234.8554 

24. 1739 
155. 1730 
781.2983 
238.0367 

23.8494 

155.3176 

799.5397 

237.2727 

22.8404 

158.4916 

975.4359 

241.0397 


B-2 


, 685.7206 

382.2241. 

234.5487 

636.6875 

224,5776 

181.4015^ 

140.01S7  - 

213.8|84 

24.6692 

24.0958 

30.1119 

23.8058 

150,7653 

142.7408 

142.3283 

452.  W 4 

582.8910 

452.8691 

99.3572 

692.5345 

226. 084 J 

177.6558 

156.94^4 

229.6438 

24,0614 

24. 9972 

30.1761 

25.3786 

151.3428 

139.5786 

141.7217 

152.0469 

631 . 1062 

439. 1368 

112. 1004 

603.434^6 

225.9728 

154.2619 

141.1226 

233.7623 

21.2186 

24.9119 

27.7333 

23.7814 

150. 1872 

141.8947 

141.9033 

153.2708 

732. 1046 

408.3553 

242. 1310 

717.1751 

213.0835 

170.9062 

165.0574 

i 231.4337 

21.6078 

24.9439 

27.0881 

23.4747 

150.0167 

142.2869 

144.0564 

154.8847 

707 . 1 235 

403.7977 

277.7309 

795.4900 

213. 3606 

173.8675 

190.8575 

234.9182 

24. 1544 

26.8138 

^ * 27.4800 

22.9625 

148.4039 

141.3931 

145.6253 

153.7467 

530.2255 

302.4558 

282.6504 

770.1820 

213.4558 

162.6229 

208.6913 

1 230.1079 

23.2053 

27.9811 

28.0453 

1 

22.5064 

146, 5494 

i 141.3203 

146.3633 

157.9086 

538.5653 

237.7077 

275.0908 

961.7020 

200.8970 

156.9168 

219.2526 

238.9160 

25. 7386 

27.6131 

25.6611 

145.9717 

142.9464 

148.7797 

385.2754 

241.5809 

467.4952 

203.9913 

178.3382 

221.3414 

25.5519 

28.4825 

24.3919 

146.2350 

140.7200 

150.3806 

401. 1778 

225.4498 

581.5632 

205.2599 

146.5065 

1 

223.4382 

24.6292 

1 28.3597 

23. 1225 

144.7483 

1 140.7383 

150.3558 

. 

429.6817 

231.2071 

640.9430 

192. 2712 

147. 6860 

XongitudfH*  1 5 . 1 69.9 
rang*  •706.1149 
ing  »23!&.3438 

24.0845 

156.4569 

784.3776 

230.2172 

24.1191 

155.2021 

729.4729 

235.0567 

23.5569 
154.2322 
712,7239 
230. 1922 

21.7024 

151.4458 

703.6050 

214.6631 

24,3380 

151.7335 

584.6971 

224.1484 

26. 1924 
152.6832 
545. 1220 
236.0152 

25.9971 

153.0586 

568.0233 

236.3263 

25.4924 

148.5497 

426.6414 

212.7171 

26.8048 
149.8825 
407.8.640 
£27. 1970 

22.8895 

147.8464 

552.6541 

200.8624 


148, 1508 

i 143.:4494 

387.5504 

251.4282 

212.4499 

9.4871 

23.3131 

35.6381 

147.2220 

145.3606 

517.8066 

251.2328  : 

198.1862 

348.0971 

24.1918 

35.3680 

145.4252 

. 143.8789 

445.2249 

260.2054 

187.9807 

4.5052 

27.5042 

35.6175 

147.0110 

142.3557 

280.7504 

262.9130 

210.9633 

21 . 1484 

25.9006 

: 36.5841 

141.0524 

i 142.5895 

379.2604 

314.0235 

152.4379 

27.6159 

143.2209 

242.3735 

166.3115 

28.7610 

141.5502 

219.6241 

138.8161 

29.4627 
145. 1253 
132.0519 
199.0869 

31 . 6398 
144.8554 

29.0025 
281 . 4249 

27. 7574 
142.4706 
246.4452 
156.7858 

33.4342 

145.2125 

122.4009 

338.0502 


15. 2238 

37.6162 

144.3749 

364.2763 

359.4306 

37. 8918 
143.4401 
383, 1361 
6. 104? 

39.0103 

144.0451 

448.0725 

1.5173 

38,9919 

143.4910 

448.5386 

4.8320 

35.5481 
144.7190 
241 . 0876 
355. 1162 

35.3524 
147, 1819 
270. 1666 
328. 4981 


147,5337 

166,61^7 

264.7308 

31.6834 

145.9668 

,85.6266 

275.9927 

29.2026 

154.9421 

568.3192 

258.4252 

25.2766 

150.8612 

511.0254 

224.2479 

21.6056 

157.2165 

913.0576 

232.3567 

24.2772 

152.0889 

600.2352 

225.3372 

24.7600 

146.5269 

423.7997 

196.6842 

24.9467 

145.8947 

404.7838 

192.4123 

23.6586 

142.3261 

484.6538 

167.0323 

23.8719 
142.5450 
469.7064 
168. 1423 

24.9225 
142. 5S44 
407. 5494 
166. 7382 
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23.903^ 
142. -86 1 7 
464.7772 
170.2050 

24.9439 

141.7767 

417.8726 

160.7874 

25.4292 

141.2675 

400.2164 

155.7217 

25.2853 

141.4722 

404.0835 

157.6524 

28.2400 

141.8617 

235.3568 

146.7951 

28.5919 
141.9769 
214.3119 
145. 1798 

27. 1892 
144.2258 
261.3946 
179. 1445 

27.0267 

147.7375 

325.3464 

214.4299 

24.0319 

150.6431 

562.4042 

218.2910 

22.0239 
156.3961 
862. 6623 
23 i . 4802 


21.4664 
151.4708 
716.4397 
214. 1641 

20.2292 

149.6703 

737.9595 

204.3521 

19.9944 

148.0131 

721.3589 

196.9925 

18.6106 

144.9367 

776.8514 

182.6986 

20.8400 

142.7089 

647.9674 

172.0449 

23.3786 

143.4169 

492.2376 

174.3178 

23.4986 
142.9447 
488.6098 
171 /1987 


23. 1067 
143. 1406 
510.0554 
172.7742 

22.3439 

142.8181 

557.7144 

171.4900 


22.7492 

550.21^7 
197. 1376 

2,1.7706 

149.4092 

647.0060 

206.2359 

21.7172 

150.0683 

665.6214 

209.0362 
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FIGURE  Cl.  (U)  Number  of 
Co’^ rtsi  (upper  number)  r 
Pr»‘.p  x t l''  1 of  Area  Covered 
(m:^  I imber)  on  10  September 
for  X 5*  square.  Ship 
Density  is  not  calculated  if 
less  than  half  the  square  was 
covered.  (U) 
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FIGURE  C2.  (U)  Number  of  Contacts 

(upp^  nuxober) , Proportion  of 
Area  Covered  (middle  number)  on 
12  September  for  5**  x 5”  square. 
Ship  Density  is  not  calculated 
if  less  than  half  the  square  was 
covered.  (U) 
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FIGURE  C3.  (U)  Number  of 
Contacts  (upper  number) , 
Proportion  of  Area  Covered 
(middle  number)  and  Ship 
Density  (lower  number)  on 
14  September  for  5*  x 5* 
square.  Ship  Density  is 
not  calculated  if  less  than 
half  the  square  was  covered. 
(U) 
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FIGURE  C4.  (U)  Ntunber  of 
Contacts  (upper  number) , 
Proportion  of  Area  Covered 
(middle  number)  and  Ship 
Density  (lower  nximber)  on 
16  Septen^er  for  5®  x 5® 
square.  Ship  Density  is 
not  calculated  if  less  than 
half  the  square  was  covered. 
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figure  C6.  (U)  Nuiaber  of  Contact 
(upper  number) , Proportion  of 
Area  Covered  (middle  number) 
and  Ship  Density  (lower  numbez 
on  20  September  for  5*  x 5* 
square.  Ship  Density  is  not 
calculated  if  less  than  half 
the  square  was  covered.  (U) 


FIGURE  C7 . (u)  Number  of  Contacts] 
(upper  number) , Proportion  of  j 
Area  Covered  (middle  number)  and 
Ship  Density  (lower  number)  on 
22  SeptOBbei  for  5*  x 5*^  square. 
Ship  Density  is  not  calculated 
if  less  than  half  the  square 
was  covered.  (U)  j 
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APPENDIX  B.l 


SHIP  DENSITY  DETERMINATION  WITH 
HIGH  RESOLUTION  SKYWAVB  RADAR  SURVEILLANCE: 
REPORT  OF  FEASIBILITY  EXPERIMENT  (U) 


1.  (C)  INTRODUCTION  (U) 

(C)  The  Hide  Aperture  Research  Facility  (WARF)  hlgh~resolution 

OTH-B  radar  was  operated  over  the  period  of  10  Sept^ber  1975 
through  14  September  1975  for  the  purpose  of  detecting  ships  in 
the  Pacific  Ocean.  This  mission  was  Identical  to  that  described 
in  the  CHURCH  OPAL  Exercise  Plan.  The  total  area  of  coverage  is 
illustrated  in  Figure  E.1-1.  Daily  ship  surveillance  was  per- 
formed in  areas  exceeding  5**  latitude  x 5*  longitude  that  were 
centered  on  the  following  locations: 

10  September  1975;  35*  - 40*N  and  140*  - 145*W 

11,12  September  1975:  30*  - 35*N  and  140*  - 145*W 

13,14  September  1975:  25*  - 30*N  and  140*  - 145*W 

The  radar  was  operated  for  a period  to  time  close  to  10  hours  on 
each  day. 

2.  (C)  RESULTS  (U) 

2.1  (C)  DATA  SUMMARY  (U) 

(C)  Approximately  55,000  individual  ocean  surface  cells  were 

sampled  each  day  for  the  purpose  of  detecting  ships.  Each  cell 

measured  7.5  km  in  radar  range  depth  and  approximately  13  km  in 

2 

radar  azimuth  width,  or  about  100  (km)  per  cell.  The  ocean  area 

2 

surveyed  per  day  measured  about  400,000  (km)  . This  area  accounts 
for  4,000  of  the  55,000  sampled  cells.  On  the  average,  over  half 
of  the  radar  dwells  were  usable  on  each  day;  thus,  each  5*  x 5* 
area  was  seunpled  several  times  (some  portions  more  than  others) . 
When  necessary,  many  dwells  were  devoted  to  the  verification  of 
ship  targets  by  concentrated  sampling  of  single  areas.  The  un- 
usable radar  dwells  contained  meteor  echoes  (which  camouflage 
ships) , insufficient  signal  strength  (necessitating  radar  fre- 
quency changes)  or  unusually  spread  sea  clutter  due  to  disturbed 

E-1 
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ionosphwic  conditicms. 

(S)  Suimatirles  of  the  bcefln  areAs  aurveyedi  on  each  day  with 

good  data  are  presented  in  Figures  E.1-2  through  E.I*>6.  'fhe  area 

corresponding  to  each  radar  dwell  is  outlined  (approximately 

rectangular) . Each  of  these  dwells  contained  105  individual  (but 

2 

slightly  overlapping)  100  (Ion)  cells  that  were  recorded  simul^ 
taneously . Some  areas  were  sampled  more,  often  and  thus  have  higher 
probability  for  determination  of  ship  density.  (The  number  of 
overlapping  dwells  did  not  reproduce  in  the  figures.)  Approxi- 
mately 4 good  dwells  on  an  area  were  necessary  in  this  experiment 
to  insure  detection  of  ships  of  approximately  400  feet  and  larger, 
at  radial  speeds  in  excess  of  9 or  10  knots.  Smaller  ships  may  or 
may  not  have  been  detected. 

2.2  (S)  SHIP  TRACKS  (U) 

(S)  The  tracks  of  ships  detected  on  the  5 days  are  presented 

in  Figures  E.1-7  to  E.1-37,  in  chronological  order  measured  to 
the  beginning  of  each  track.  The  date  of  the  track  is  Printed  at 
the  top  of  each  map.  (To  simplify  software  developed  recently 
for  this  test  the  maps  were  produced  with  varying  scale  ratios  for 
latitude  and  longitude.)  Except  where  noted,  each  ship  track  is 
drawn  on  the  map  with  an  arrow  denoting  the  direction  of  travel. 
Data  pertinent  to  each  track  are  printed  on  each  figure. 

• The  top  line  gives  the  begin  and  end  time  (GMT) 
of  the  track,  ihe  track  duration,  and  the  begin 
and  end  positions  in  decimal  degrees) . 

• The  third  line  gives  the  average  observed  radial 
speed  (kts) . 

• The  fourth  line  gives  the  average  estimated  radar 
cross  section  (AVG  RCS) . 

• The  fifth  line  gives  the  calculated  true  speed  (kts) 
and  course  heading  at  the  end  of  the  track  * 

For  a few  tracks  (e.g..  Figure  E.1-9)  the  true  speed  and  heading 
has  been  omitted  because  either  too  few  hits  were  present  or  the 


time  span  too  short  to  establish  an  accurate  cross*ratnge  rate, 
in  those  cases  we  have  Xat>eled  each  track  as  "Inbound"  or  "out-* 
bound" , relative  to  central  California  t based  on  the  sign  of  the 
radial  speed.  The  accuracy  of  the  track  end  points  is  within 
20  km  (10  nml) . 

2.2.1  (s)  Probability  Estimates  (U) 

(S)  A probability  has  been  associated  with  each  track  and 

is  printed  on  the  lower  right-hand  corner  of  each  figure.  This 
probability  is  based  on  the  number,  time  spacing,  and  consistency 
of  radar  contacts.  These  probabilities  are  rough  estimates,  and 
may  be  defined  as  follows: 

e >90%  means  the  track  is  almost  undoubtedly  a 
ship. 

e >50%  means  that  something  less  than  half  of 
all  the  tracks  with  this  probability  are 
probably  only  coincidental  correlations  of 
false  echoes  produced  by  noise,  meteors, 
clutter  enhancements,  and/or  aircraft. 

• <25%  means  that  probably  less  than  one  fourth 
of  the  tracks  with  this  probability  are  ships. 

In  some  hypothetical  strategic  situations,  the  radar  operator 
would  surely  pass  all  the  90%  tracks,  and  probably  the  50%  tracks. 
Further  sensor  verification  of  the  25%  tracks  would  be  necessary 
to  give  them  reasonable  validity. 

2.2.2  (S)  Track  Combinations  and  Multiple  Targets  (U) 

(C)  In  two  cases.  Figures  E.1-21  and  E.1-31  (each  in  three 

parts),  we  had  two  tracks  spaced  by  two  and  seven  hours,  respec- 
tively, that  could  be  combined  into  one  track.  In  each  case,  this 
conclusion  was  based  on  the  consistency  of  both  the  target's  radial 
speed  and  the  final  speed  estimate. 

(U)  For  Figure  E.1-21  we  would  use  only  the  combination 

track  (Figure  E.1-21-C). 
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(S)  Figure  E.1-31  is  much  more  complicated.  We  have  high 

probability  for  both  the  begin  and  end  tracks  (spaced  by  seven 
hours) « and  very  good  consistency  for  the  final  track.  However, 
we  probably  had  more  than  one  ship  in  a group  of  2 or  3 (spaced 
by  tens  of  km)  in  each  of  the  tvro  tracks.  Whether  or  not  the  sane 
ships  that  were  together  on  the  first  track  (Figure  E.1'31“A) 
were  also  together  on  the  second  track  (Figure  E.1-31-B)  is  uncer*' 
tain.  Owing  to  the  long  time  between  tracks,  we  would  recommend 
that  the  two  tracks  be  left  separate,  with  two  ships  per  track. 
Figure  E.1-31-C  would  be  ignored  unless  other  sensor  information 
can  help  relate  the  two  separate  tracks. 

(C)  The  detection  log  that  formed  the  track  of  Figure  E.1-8 

also  suggested  the  presence  of  two  or  three  closely  spaced  ships. 

We  would  recommend  the  association  of  two  ships  for  this  track. 

2.2.3  (S)  Target  Size  (U) 

(S)  The  ship  size  is  roughly  related  to  its  average  RCS 

estimated  from  the  data  and  printed  beneath  each  track.  This  esti- 
mate is  based  on  an  estimate  of  the  sea  clutter  cross  section  which 
has  been  found  to  be  surprisingly  predictable.  Based  on  measure- 
ments of  both  full-size  and  scale  model  ship  cross  sections,  we 
would  divide  the  ship  sizes  as  follows: 

• Small  (<400  feet)  for  RCS<35 

9 £4edium  (400  to  600  feet)  for  35£RCS^40  dB  ^ 

2 ® 

e Large  (>600  feet)  for  RCS>40  dB^ 

3.  (S)  CONCLUSIONS  AND  RECOMMENDATIONS  (1) 

(S)  A total  Of  32  tracks  (containing  a total  of  34  ships  were 

derived  by  means  of  surveillance  by  high  resolution  OTH-B  radar 
during  the  period  of  10-14  September  1975.  The  surveillance  areas 
are  outlined  in  Section  1 of  this  report.  Of  these  tracks  we 
attach  very  high  confidence  (>90%)  and  medium  confidence  (>50%) 
as  follows: 
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SECRET 


''ll 


(S) 


10  Sept«i>ber 


(>90%) 

(>50%) 


TF* 


TOTAL 


2 tracks,  3 ships 
2 tracks,  2 ships 

4 tracks,  5 ships 


11  September 


2 tracks,  2 ships  (>90%) 
2 tracks,  2 ships  (>50%) 


TOTAL  4 tracks,  4 ships 


12  September 


TOTAL 


1 track,  1 ship  (>90%) 

1 track,  1 ship  (>50%) 

2 tracks,  2 ships 


13  September 


4 tracks,  4 ships  (>90%) 
2 tracks,  2 ships  (>50%) 


TOTAL 

14  September 


6 tracks,  6 ships 

5 tracks,  7 ships  (>90%) 
2 tracks,  2 ships  (>S0%) 


TOTAL  7 tracks,  9 ships 


On  10  September  one  track  had  two  ships  closely  spaced,  and  on 
14  September  two  tracks  had  two  ships  in  a group.  Tracks  remaining 
of  the  34  had  a confidence  of  less  than  25%  and  are  included  herein 
only  for  academic  purposes  (e.g.,  if  another  sensor  could  verify 
them) . The  position  accuracies  for  all  tracks  are  20  km  (10  nmi) . 

(S)  Ship  radial  speeds  ranged  from  10.8  to  22.6  kts,  while 

estimated  true  ship  speeds  ranged  from  14.0  to  24.6  kts.  Both 
Inbound  and  outbound  ships  were  tracked.  The  accuracy  of  the  ra- 
dial speed  is  about  5%.  The  course  and  true  speed  accuracy  is  pri- 
marily determined  by  the  time  span  between  track  end  points  when 
the  20  km  position  accuracy  is  applied. 

(C)  Several  new  area-surveillance  ship  tracking  methods  were 

developed  as  a result  of  this  feasibility  experiment.  Heret  ^ore, 
(with  one  exception)  we  have  employed  WARF  for  the  tracking  only 
cf  specific  (single)  ships  of  high  interest  either  to  SRi  or  the 


U.S.  Fleet.  The  surveillance  of  relatively  large  areas  with  high 
resolution  i quires  rather  careful  planning  to  insure  efficient 
target  hit-to-hit  correlation.  New  correlation  procedures  were 
developed  in  the  CHURCH  OPAL  exercise,  and  better  software  is 
now  being  written  to  autcxnate  much  of  what  was  done  manually  after 
the  experiment  (this  work  is  part  of  a separate  ONR-supported  task) . 
Additionally,  %re  are  now  prepared  to  suggest  new  scanning  proce- 
dures, possibly  a fence  concept,  for  any  future  experiments  with 
LIU^P.  A suitable  blend  of  scanning  strategy  and  hit -to-hit  corre- 
lation procedures  %#ould  yield  target  tracks  all  of  high  confidence. 
Tracks  of  initially  lower  confidence  vx>uld  be  verified. 
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FILES  11^  151  TIME  16’41«O0  19-.£i6*e0  3.42  HOURS  28.74  142.75 
SPEED  RANGE  13.8  15.5  HIM  SNR  13  0 28.93  142.12 
AUG  SPEEO=  14.1  SPEED  FROM  LSQ^  10.4 

AUG  SNR«  18.5  AUG  RCS«  35.6  RCS  S.D.M*  1.4  AUG  SPEED 
SPEED  FROM  LAT/LOH^  14.5  HEADING  71 . 1 DEGREES 


>90% 


Figure  E.1-22 
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31.601 


30.56 


30.09 


FILES  117  252  TIME  18:32  = 08  21:33*41 
SPEED  RANGE  13.9  14.3  HIH  14.2 
AUG  SPEED-  14.1  SPEED  KROtl  LSQ»  14.2 
AUG  SHR»  14.7  AUG  RCS-  32.7  RCS  8.0.M«  -.3 
SPEED  FROM  LAT/LON-  26.8  HEADING  198.7  DECREES 


3.13  HOURS  29.55  142.63 
29.20  141.46 


AUG  SPEED 


< 25% 


Figure  E.1-24 
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144.00  143.25  142.50  141.75  141.00 

FILES  191  513  TIME  19 « 49 « 00  01«45«45  5.93  HOURS  29.43  141 .44 

SPEED  RANGE  -14.5  -11.7  MIH  SHR  11.8  27.76  142.00 

AUG  SPEEOa.  -12.4  SPEED  FROM  LSQ«  -13  6 

AUG  SNR«  13.5  AUG  RCS«  43.2  RCS  S.D.M=:  -1.6  AUG  SPEED 
SPEED  FROM  LAT^LON-  17.7  HEADING  196.4  DEGREES 


HI 


> 90% 


Figure  E.1-25 
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27.50 


27.00 


26.50 
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147.00  146.25  145.50  144.75  144.00 

FILES  302  556  TIME  22  30 ‘24  02«27‘05  3.94  HOURS  26.56  145.34 

SPEED  RANGE  15.3  15.8  MIN  SNR  13.6  28.00  144.94 

AUG  SPEED*  15.5  SPEED  FROM  LSQ*  24.4 

AUG  SNR«  15.9  AUG  RCS*  32.1  RCS  S.D.M*  -.7  AUG  SPEED 
SPEED  FROM  LAT/LON»  22.7  HEADING  13.7  DEGREES 


Figure  E.1-26 


> 60% 


E-34 


UNCLASSIFIED 


SEP  13  1975 
29.o0r 


28.50 


28.68 


27.50 


27.001 

145.00 


UNCLASSIFIED 


« <1HB0UMD) 


144.50 


144.00 
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FILES  387  488  TIME  23<48<27  01>23>59  1.59  HOURS  28.26  144.43 

SPEED  RANGE  15.6  16.0  NIN  SHR  13.4  27.96  143.70 

* AUG  SPEED-  15.8  SPEED  FROM  LSQ-  10.2 
AUG  SNR»  13.9  AUG  RCS-  33.5  RC8  S.D.M-  -2.6  AUG  SPEED 


< 351 


Figure  E.1-27 
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142.50 


142.25 
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Files  425  462  time  ee»25«e0  e0<59>42  .58  hours  29.67  142  32 
SPEED  RANGE  18,0  11,7  MIN  SNR  1-3.9  29.57  142  l£ 
AviC  SPEEO=  11.1  SPEED  FROM  LSQ»  14.2 

j5^  SNg®  16.2  AUG  RCS*  38.7  RCS  S.D.Ma  -2.0  AUG  SPEED 
SPEED  FROM  LAT/LON=  21.2  HEfiDING  120.4  DECREES 


> 80% 


Figure  E.1-28 
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144.09 
FILES  7 


143.75  143.50  143.23 

38  TIME  16:3S>’00  I7»09*3l 

SPEED  RANGE  9.7  12.7  MIH  SHR  13.0 

• AUG  SPEED-  11.0  SPEED  FROM  LSQ*  .6 
AUG  SNR«  15.6  AUG  RCS=  40  3 RCS  S.D.M- 


143.00 

.51  HOURS  30.22  143.79 
30.41  143.76 


-1.1  AUG  SPEED 


1 to  3 ships  in  group 
> 90% 


Figure  E.1-31-A 
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FILES  422  482  TIME  2'i-47’ed  ee>48*22  1 .82  HOURS  38.58  141 .93 
SPEED  RANGE  9.7  12.4  MIN  SMR  13.6  30.48  141.65 
AUG  SPEED**  16.8  SPEED  FROM  LSQ-  21.4 

AUG  SNRa  19.5  * AUG  RCS**  36.8  RCS  S.D.M»  -.9  AUG  SPEED 
SPEED  FROM  LAT/LON»  15.5  HEADING  189.6  DEGREES 


1 or  2 ships  in  group  > 90% 


Figure  E.1-31-B 


(sea  7-482) 
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144.00  143.25  142.50  141.75  141.00 

FILES  7 482  TIME  16^39*00  00*48*22  8.16  HOURS  30.39  143.65 

SPEED  RANGE  9.7  12  7 Mill  SHR  13.7  30.40  141.76 

AUG  SPEED**  10.8  SPEED  FROI1  LSQ*=  11.2 

AUG  SNR»  18.5  ‘ AUG  RCS»  38.3  RCS  S.D.M*  .0  AUG  SPEED 
SPEED  FROM  LAT/LON-  12.8  HEADING  5>O.0  OEGP.EES 


1 to  3 ships  in  group 


>90% 


(combination  of  last  2 tra 


Figure  E.1-31-C 
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FILES  62  498  TIME  17:31 « 25  01*02*21  7.52  HOURS  26.18  140.20 
SPEED  RANGE  10.4  13.6  MIH  SI{R  13.3  28.68  140.10 
AUG  SPEED*  11.5  SPEED  FROM  LSQ*  10.9 

AUG  SNR*  13.4  AUG  RCS*  35.3  RCS  S.D.M=  -2.1  AUG  SPEED 
SPEED  FROM  LAT/LOH*  28.3  HEADING  2.0  DEGREES 


Figure  E.1-32 


E-42 


UNCLASSIFIED 


> 50% 


UIWLAIiSrNtU 


H SEPT  X9?5 

ji.eeP 


30.25 


29.50 


28.75 


28.80L- 

i47.ee 


146.50 


14b. 08 


145.50 


145.00 


FILES  106  407  TIME  18*19524  23-29*50  3.17  HOI 

SPEED  RANGE  -17.3  -14.0  MIN  SNR  12.3 
AVG  SPEED*  -15.5  SPEED  FROM  LSQ*  -12.0 
AUG  SNRa  15.5  AUG  RCS=  34.9  RCS  S.D.M*  1.7 
SPEED  FROM  LAT/LON*  18.9  HEADING  279.0  DEGREES 


3.17  HOURS  28.87  145.09 
29.11  146.93 


AUG  SPEED 
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146.50 


145.25  145.00 


FILES  316  399  TIME  22-06:00  23«22*53  1,28  1 

SPEED  RANGE  11.9  13.8  MIN  SNR  12.0 

AUG  SPEEDa  12.9  SPEED  FROM  LSQa  15.2 
AUG  SNR-  13.0  AUG  RCS-  24.5  RCS  S.O.M- 
SPEED  FROM  LAT/LONa  14.0  HEADING  87.1  DEGREES 


1,28  HOURS  29.19  146.55 
29.21  146.21 


-.5  AUG  SPEED 


> 50% 


Figure  E.1-34 
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FILES  321  387  TIME  22- U -09  23>12<38  1 .02  HOURS  27.76  145.81 
SPEED  RANGE  -17. S -16.8  MIH  SHR  12.6  27.66  145.33 
AUG  SPEED*  -17,1  SPEED  FROM  LSQ*  -8  2 

AUG  SNR-  18.5  AUG  RCS*  32.4  RCS  S.D.M-  -.6  AUG  SPEED 
SPEED  FROM  LAT/LOM*  17.6  HEADING  251.7  DEGREES 


Figure  E.1-35 
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FILES  322  371  TIME  22^12^00  22»58«22  .77  HOURS  26.70  144.97 
SPEED  RANGE  19.2  20.7  MIN  SJIR  15.0  26.82  144.65 
AUG  SPEED*  20.1  SPEED  FROM  LSQ*  23.8 

AUG  SNR*  25.6  AUG  RCS*  34.2  RHS  F.D.M*  .7  AUG  SPEED 
SPEED  FROM  LAT^LON*  22.2  HEADING  07  2 DEGREES 


Figure  E.1-36 
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APPENDIX  E.2 

OTH  RADAR  SHIP-DISTRIBUTION  MEASUREMENTS:  SEPTEMBER  1975 

1.  (U)  INTRODUCTION  (U) 

(U)  Surface  shipping  surveillance  tests  were  conducted  with 

the  SEA  ECHO  over-the-horizon  (OTH)  zedar  at  San  Clement  Island 
during  the  period  from  September  20-24,  1975,  for  a project  spon- ' 
sored  by  the  Office  of  Naval  Research.  The  survey  was  made  by 
daily  scans  of  specified  areas  In  the  Gulf  of  Alaska  with  the  ob- 
jective of  defining  shipping  density  distributions  in  support  of 
analyses  being  conducted  by  Planning  Systems  Incorporated  of 
McLean,  Virginia.  This  report  gives  the  results  of  this  investi- 
gation together  with  a description  and  some  discussion  of  the  ex- 
perimental parameters. 

2.  (C)  TASK  OBJECTIVES  (U) 

(C)  The  objective  of  the  OTH  surveillance  tests  was  to  pro- 

vide information  on  the  density  distribution  of  surface  shipping 
in  a specified  area  and  for  specified  times  for  comparison  and 
use  in  theoretical  models  for  ship  distribution.  More  specifi- 
cally the  scope  of  this  specific  work  was  categorized  as  a test 
of  the  capability  of  OTH  radar  systems  to  provide  ship  surveil- 
lance data  of  adequate  utility  for  such  comparison.  The  test  re- 
quirements called  for  a dally  snapshot  of  specified  5*  x 5*  (lat.- 
long.)  areas  to  be  performed  with  a spatial  resolution  which  should 
be  adequate  for  a reasonable  detail  of  target  distribution  within 
this  area,  and  an  exposure  time  was  to  be  consistent  with  a rea- 
sonable freezing  of  15  kt  ship  motion.  These  objectives  were 
achieved . 

3.  (C)  TEST  DETAILS  (U) 

(C)  Two  test  areas  were  specified,  namely  a region  bounded  by 

140**  and  145**W  and  by  45**  and  50**N,  and  a second  region  bounded  by 
140**  and  145**W  and  by  50**  and  55**N.  The  first  of  these.  Area  1, 
was  specified  for  surveillance  on  September  22  and  23,  and  the 
second.  Area  2,  on  September  20,  21  and  24.  By  agreement  with  PSI, 


iry- 


the  dally  ti^  of  eaoh  gurveillanca  — for  ,thia  initial  exercise  — 
was  to  be  determined  at  the  SEA  ECHO  site  on  a or^erion  of  best-- 
ionospheric  propagation  oondltions.  Accordingly  the  surveiliance 
scans  were  conducted  generally  at  mid*^ay,  Data  were  taHen  for 
the  specified  areas  on  Septonber  21,  22  and  24. 

(U)  Systoa  difficulties  prohibited  usg^ful  sieasureiiients  on 

September  20  and  23.  ^hfo  scans  are  ^eiienti^  for ^'dif iterant  times, 
towevrer,  on  Septeodder  21. 

(U)'  The  SEA  ECHO  radar  was  set  up  for  the  maximum  spatial 

resolution  (smallest  cell  size)  for  maximum  signal-to-clutter  ratio. 
The  ;ell  size  used  was  a nominal  2**  in  beam  width  and  15  km  (8.09  nm) 
in  range.  The  data  is  reported  on  the  basis  of  this  resolution. 

(Note : Zn  the  accompanying  graphical  representations,  for  con- 
venience in  plotting  the  range-wise  grid  size  is  shown  in  steps 
of  10  nm.)  In  presenting  the  data,  each  target  is  shown  in  the 
center  of  each  resolution  cell.  However,  no  further  resolution 
can  be  assigned  than  the  presence  of  a target  in  a cell.  (Identi- 
fication of  more  than  one  target  per  cell  is  readily  possible  if 
their  doppler  signatures  differ,  however.) 

(U)  In  a scan  of  a given  area,  20  consecutive  range  cells  were 

interrogated  simultaneously  in  successive  azimuthal  beam  positions. 
The  system  was  then  reset  for  the  second  (more  distant)  set  of  20 
range  cells  and  again  sequenced  at  these  ranges  through  the  suc- 
cessive six  beams.  Each  20-cell,  six-beam  scan  re<]uired  about  10 
minutes.  Re-set  to  the  new  ranges  also  required  about  10  minutes. 
Therefore  an  entire  scan  was  accomplished  in  30  minutes. 

(C)  Since  the  radar  spatial  resolution  cell  size  is  much  lar- 

ger than  that  of  Individual  ships  (both  in  physical  and  radar 
cross-section) , doppler  - or  velocity  - processing  must  be  invoked 
to  discriminate  the  targets  of  interest  from  the  sea  return  (clutter) 
For  these  tests,  the  SEA  ECHO  processor  was  arranged  to  acquire 
at  each  beam  position  for  approximately  25  seconds  producing  a 256 
point  Fourier  transform  for  each  resolution  cell  in  that  beam  be- 
fore sequencing  to  the  next  beam.  The  sum  of  four  such  transforms, 
taken  at  intervals  of  the  beam  revisit  time,  were  averaged  and 
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presented  on  a CRT  for  display.  Bach  display  presented  spectral 
plots  of  radar  cros8*>scK!tion  vs.  doppler  frequency  (target  velocity) 
for  ten  consecutive  range  celiis.  'Bach  display  WaS  photographed » 
and  the  resulting  Polaroid  print  was  used  in  a visual  examination 
for  the  presence  of  ship  targets  Ml  data  acquired  in  these 
tests  were  digitally  recorded  in  a 1.  oat  which  is  suitable  for 
more  sophisticated  machine  processing  for  discrimination  and  detec~ 
tion  of  targets.  While  such  pxoceosing  is  beyond  the  scope  of 
the  tasks  defined  in  this  preliminary  work,  it  is  recommended. 
Details  of  this  and  the  data  par«uneters  are  discussed  in  a later 
section  of  this  report. 

4.  (8)  RESULTS  (U) 

(U)  Results  are  displayed  in  the  attached  charts  each  of  which 

is  identified  by  area,  date  and  time.  These  charts  are  to  scale. 
Latitude  and  longitude  in  truee  coordinates  are  given  in  a 1 x 1* 
grid.  Radar  coordinates  are  given  in  nautical  miles  and  true  bear- 
ings from  the  SEA  ECHO  site.  The  radar  azimuthal  grid  is  given 
in  a 2*  resolution  corresponding  to  the  half -power  points  of  the 
radiated  beams  — i.e.,  each  beam  lies  within  adjacent  labeled 
grid  lines. 

(S)  Bach  target  identified  is  shown  by  an  arrow  whose  direc- 

tion indicates  whether  the  velocity  vector  was  generally  away  from 
or  toward  the  radar.  Further,  a discrimination  is  shown  between 

targets  whose  radar  cross  section  is  estimated  to  be  either  greater 

2 

or  less  than  10, DOOM  , i.e.,  40  dB  greater  than  on  square  meter 
(40  dBsm) . A double-line  arrow  indicates  a radar  cross  section 
which  is  greater  than  40  dBsm  and  one  with  a single  line,  a tar- 
get smaller  than  40  dBsm.  (For  reference,  the  cross-section  of  a 
large  U.S.  Navy  destroyer  is  approximately  33  dBsm  at  HF  wave 
lengths . ) 

5.  (S)  DISCUSSION  (U) 

(C)  A large  ship  (i.e.,  larger  than  a destroyer  but  smaller 

than  an  aircraft  carrier)  can  be  generally  categorized  in  the 
10,000  square  meter  HF  radar  cross  section  range. 


(S)  The  competing  signal  return  from  the  minimum  ra^dai^  resolu’ 

tlon  cell  (at  a range  of  1500  nm)  Is  of  the  order  . of  1.4  x 10 V. 
Doppler  processing  Is  utilized  to  extract  targets  of  interest  from 
this  significantly  larger  clutter.  The  following  details  are 
pertinent . 

(S)  The  spectra  produced  by  the  doppler  processing  represent 

a continuous  set  of  values  of  radar  cross-section  «S,  a function  of 
doppler  frequency  (target  velocity) . The  SEA  ECHO  display  has 
been  Implemented  to  show  both  receding  and  approaching  doppler . 
(Thus  processed/  clutter  return  from  the  sea  always  appears  as  two 
spectral  lines  symetrlcally  disposed  about  zero  doppler  frequency. 
These  lines  represent  sea  waves  receding  and  approaching  the  radar 
Because  of  spatial  ("Bragg")  resonance/  the  waves  producing  these 
spectral  lines  are  those  components  whose  wave  length  Is  exactly 
half  of  the  electromagnetic  wave  length  corresponding  to  the  fre- 
quency selected  for  the  radar  operation.)  Furthermore,  because  of 
the  dispersion  relation  between  ocean  wave  lengths  and  their 
feloclty,  the  doppler  frequency  (velocity)  produced  by  these  waves 
are  exactly  determined.  For  example,  at  the  (approximate)  SEA 
ECHO  frequency  of  16.0  MHZ  used  in  these  tests  the  sea  resonant 
spectral  (Braff)  lines  peak  at  ±0.41  Hz,  representing  a radial 
velocity  of  ±7.44  kts.  Hard  targets  such  as  ships,  whose  radial 
velocity  is  at  or  near  this  value,  will  bo  masked  by  the  effects 
of  clutter,  particularly  in  the  case  of  the  dominant  Braff  line, 
i.e.,  the  one  which  is  being  produced  by  the  wind -driven  sea 
which  produces  resultingly  high  cross-section. 

(C)  Choice  of  operation  with  two  or  more  simultaneous  radar 

frequencies  (a  process  which  Is  uniquely  possible  with  the  SEA 
ECHO  system)  will  produce  target  spectra  in  which  the  Bragg  lines 
are  shifted  relative  to  the  doppler  produced  by  a hard  target. 

Use  of  this  technique,  while  straightforward,  requires  some  opera- 
tionally subtle  strategies  not  invoked  in  the  limited  scope  of 
these  initial  tests.  (Braff  line  widths  are  determined  in  part 
by  the  Fourier  data  integration  period,  and  this  in  turn  must  be 
chosen  with  care  lest  doppler  smearing  is  Introduced  by  movement 
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of  or  multi-j>at'h  effects  in  thelonosphere.)  In  short,  in  the 
somewhat  simpXxfied  data  processing  used  here,  lO^M^  (and  somewhat 
smaller)  targets  are  detectable  in  the  pro^ence  off  cltitter  if 
their  doppler  signatures  (radial  velocities)  are  sufficiently 
different  fr<»a  the  clutter  doppler. 

(U)  In  the  data  presented  here,  analysis  criteria . (necessarily 

somewhat  subjective)  were  used  which  opted  for  a minimization  of 
false  alrams,  hence  toward  a reduced  probability  of  detection. 

(S)  The  following  characteristics  represent  thresholds  either 

of  which,  in  the  analysis  used,  would  have  obscured  a target. 

2 

1.  Radar  cross  section  less  than  30  dBsm  (1000  M ). 

2.  A radial  velocity  c(»nponent  of  less  than  approxi- 
mately 13  kts.  (The  doppler  spectrum  between 
the  Bragg  lines  does  provide  a window  for  low 
velocity  targets  of  something  more  than  minimal 
cross  section,  but  this  spectral  region  was  dis- 
counted in  these  data.) 

6.  (C)  TARGET  CROSS  SECTION  (U) 

(C)  Estimation  of  the  radar  cross  section  values  for  targets 

plotted  in  the  charts  was  performed  by  comparison  with  the  princi- 
pal clutter  Bragg  line.  Given  the  condition  that  at  least  one 

O 

Bragg  line  was  produced  by  a fully-arisen  sea,  the  albedo  (c  in 
radar  parlance)  for  those  waves  has  been  shown  to  be  a fixed  value 
of  -17  dB.  Since  the  resolution  cell  area  is  known,  the  sea  thus 
provides  a built-in  calibration.  The  requirement  for  the  fully 
arisen  sea  applies  only  to  the  wave  lengths  selected  by  the  radar 
frequency.  That  is,  for  the  September  20-24  period,  full  sea 
saturation  would  have  called  for  a wind  component  of  7.4  kts  in 
a radial  direction  for  a few  hours.  To  the  degree  that  this  was 
not  the  case,  the  (absolute)  value  for  a**  would  be  greater  than 
17  dB,  and  target  cross-section  should  be  estimated  as  proportion- 
ately higher. 
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7.  (C)  R£COHMENDATIOHS  (U)  | 

(C)  Based  on  the  results  achieved  in  the  tests  described 

In  the  foregoing^  NRL  recommends  two  further  approaches  to  the  OMR 
ship-density  surveillance  task. 

1.  Re-examination  of  the  September  21-24  data  using 
machine  processing. 

All  data  %fere  recorded  in  such  digital  format  ’] 

that  complete  flexibility  in  reprocessing  the  ^ 

data  can  be  accomplished  — i.e.,  in  terms  of 
Fourier  integration  lengths  and  averaging. 

Since  a characteristic  of  ionospheric  propaga- 
tion is  a time-dependent  Faraday  rotation  of 
the  EM  plane  of  polarization,  fading  of  teurget 
signals  arises  as  a result  of  target  cross- 
section  changes  with  polarization.  (This  is 
particularly  true  for  small  craft  whose  mast- 
head height  is  critical  to  cross  section.) 

The  four-integral  100  sec.  average  time  selected 
a priori  for  the  tests  is  somewhat  in  excess  of 
the  facing  period.  It  is  proposed  that  the 
individual  25  sec.  spectra  be  separately  dis- 
played and  exeunined  for  target  appearance  on  a 

i 

shorter-time  but  enhanced  signature  basis. 

Further,  the  25-second  integration  period  was 
selected  on  a synoptic  understanding  of  ionos- 
pheric time-dependence.  It  is  proposed  that 
the  data  be  reprocessed  for  a higher  doppler 
resolution  with  the  objective  of  increasing  the 
velocity  range  of  detectable  targets. 

On  the  assumption  that  for  the  areas  surveyed 

there  are  targets  of  30  dBsm  or  greater  cross- 

section  broadly  distributed  in  velocity  from  5 

to  25  kta  along  a general  great  circle  trans-  i 

Pacific  route  to  the  Far  East,  the  refinement 

of  the  above  analysis  techniques  for  September 
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^ 1975  could  increase  detected  targets  by  an 

P amount  of  the  order  of  30%. 

2.  It  is  also  proposed  that  further  surveys 
of  the  north  Pacific  shipping  densities  be 
conducted  with  the  following  optimising 
changes  in  procedure: 

With  due  recognition  that  the  surveillance 
areas  specified  for  September  1975  were 
selected  in  consonance  with  other  related 
activities,  it  should  be  noted  that  the 
resulting  radar  gound  range  of  1200  to  1700  nm 
presents  the  maximum  difficulty  to  the  SEA 
ECHO  system  (and  OTH  radars  in  general)  in 
terms  of  clean  ionospheric  propagation. 

While  our  results  show  that  ship  detection  is 
not  prohibited  at  these  ranges,  it  is 
suggested  that  considerable  benefits  could 
accrue  to  a correlation  of  theoretical 
prediction  and  direct  measurement  if  investi- 
gations were  made  at  ranges  of  the  order  of 
700  to  1200  nm.  The  reasons  for  this  are 
beyond  the  scope  of  this  report  but  generally 
have  to  do  with  the  more  or  less  unique  layers 
of  the  ionosphere  associated  with  propagation 
to  specific  ground  ranges.  SEA  ECHO  is  unique 
in  an  antenna  design  which  capitalizes  on  the 
use  of  the  lower  layers  of  the  ionosphere  — 
with  resulting  freedom  from  the  interfering 
effects  of  multiple  propagation  paths.  This 
capability  was  confirmed  in  fact  by  a test 
survey  measurement  made  on  September  24  at 
ranges  of  the  order  of  800  nm  in  which  E-layer 
propagation  produced  extremely  clean,  narrow 
clutter  spectra. 


It  is  further  proposed  that  such  measure- 
ments be  performed  using  multiple  radar 
frequencies  to  provide  redundani:  looks  Bt 
each  range  cell  with  separated  Bragg  clutter 
lines. 

^Irdly#  by  performing  area  scans  in,  a regular 
and  continued  time  sequence,  e.g.,  at  a rate 
of  one  per  hour  for  several  hours,  it  will  be 
possible  to  apply  track  correlation  tech- 
niques. By  associating  successive  target  hits 
in  progressively  adjacent  range  and/or  azimuthal 
cells  with  the  intrinsic  doppler  velocity  in- 
formation, target  absolute  velocity  vectors  can 
be  established.  Furthermore,  target  track 
techniques  can  greatly  improve  the  ratio  of 
probability  of  detection  to  probability  of  false 
alarm. 
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DEPARTMENT  OF  THE  NAVY 
OFFICE  OF  NAVAL  RESEARCH 

800  NORTH  QUINCY  STREET  - 

ARLINGTON.  VA  22217-5660"  IN  REPLY  REFER  TO 

5510/1 
Ser  43/885 
03  Dec  03 

MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj : DECLASSIFICATION  OF  CHURCH  OPAL  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  Partial  List  of  CHURCH  OPAL  Documents 

1 . In  accordance  with  reference  (a),  a declassification  review  has  been  conducted  on  a 
number  of  classified  CHURCH  OPAL  documents. 

2.  The  CHURCH  OPAL  documents  listed  in  Part.l  of  enclosure  (1)  have  been 
downgraded  to  UNCLASSIFIED  and  have  beeji  approved  for  public  release.  These 
documents  should  be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N774)  letter  N774D/3U630173,  1 1 September  2003. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 

unlimited. 


3.  If  other  CHURCH  OPAL  documents  are  located  in  your  repositories,  their  markings 
should  be  changed  and  a copy  of  the  title  page  and  a notation  of  how  many  pages  the 
documents  contained  should  be  provided  to  Chief  of  Naval  Research  (ONR  43)  800  N. 
Quincy  Street,  Arlington,  VA  22217-5660.  This  will  enable  me  to  maintain  a master  list 
of  downgraded/declassified  CHURCH  OPAL  reports. 

4.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


PEGGY  LAMBERT 
By  direction 
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